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We determine the dynamical structure factor of the two-leg spin-| Heisenberg ladder at low 
temperatures in the regime of strong rung coupling. The dominant feature at zero temperature is 
the coherent triplon mode. We show that the lineshape of this mode broadens in a non-symmetric 
way at finite temperatures and that the degree of asymmetry increases with temperature. We also 
show that at low frequencies a temperature induced resonance akin to the Villain mode in the spin-i 
Heisenberg Ising chain emerges. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Gb 



I. INTRODUCTION 

The zero temperature behaviour of two- leg spin-^ 
Heisenberg ladders is by now well understood and has 
been analyzed by a variety of theoretical methodsi"— . 
Recently, the dynamical structure factor (DSF) has 
been measured by inelastic neutron scattering experi- 
ments for the ladder compounds La4SrioCu2404i— and 
CaCu203^^ and was found to be in excellent agreement 
with theoretical predictions at T = 0. The limit of strong 
rung coupling a — J\\/ J±_ <^ 1, see Fig. [U is particu- 
larly simple. In the limit a = 0, the ground state is 
a tensor product state of rung singlets. Excitations in- 
volve breaking one of the dimers, which leads to a finite 
gap A = J_L. A small but finite J|| gives a dispersion 
to these excitations, which are commonly referred to as 
either "magnons" or "triplons"— . We will follow the lat- 
ter terminology in this work. The triplon bandwidth is 
small compared to their gap. The dominant feature of 

S=l/2 



FIG. 1. Exchange couplings for a spin-ladder system, 
strong rung coupling limit J« ^ J± . 
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the DSF at zero temperature is a delta-function follow- 
ing the triplon dispersion. At higher energies there are 
multi-triplon continua, which for small a are weak. These 
features have been analyzed in detail in the literature^"— . 
Much less is known about the finite temperature dy- 
namics of one-dimensional quantum magnets in general 
and the two-leg ladder in particular—"—. In the limit 
of large a, the DSF for the two-leg ladder model was 
studied by means of a semiclassical analysis by Damle 



and Sachde\^. They showed that at very low temper- 
atures T <C A the T — triplon delta-function peak in 
the DSF broadens and is well described by a Lorentzian 
lineshape. This behaviour was argued to be universal for 
one-dimensional gapped antiferromagnets. Very recently 
the question of how the DSF evolves as the temperature 
is increased above the semiclassical regime has been ad- 
dressed in several models by numerical^^ and analytical 
methods2^i^. It was shown that at higher temperatures, 
but still smaller than the gap, the triplon peak is broad- 
ened in a rather asymmetric fashion. In this paper we 
calculate the DSF for a spin-ladder system (Fig. [T|) at low 
temperatures. This is a quantity of experimental interest, 
probed by inelastic neutron scattering experiments^"—. 
Our calculation is restricted to the limit of weak coupling 
between the dimers, which we treat in perturbation the- 
ory to first order in a = J\i/J± for both excitation ener- 
gies and matrix elements. 

The Hamiltonian of the spin-ladder system reads 



H — Ho + Hi, 

1 L-l 
j=0 n=0 



(1) 



Here the dominant exchange coupling J± is along the 
rungs connecting neighbouring spins on different legs of 
the ladder and Jii <C J_l represents a small interaction 
between the neighbouring rungs. In the limit of zero 
interrung coupling, the ground state is a product of sin- 
glet states on every rung. The elementary excitations 
are S* = 1 triplets of energy J±, which is the difference 
between the dimer triplet and singlet states. 

Our first goal is to calculate the dynamical suscepti- 
bility, which is related to the DSF by 
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5"''(c.,Q) = — 



1 



TT 1 — exp{—/3uj) 



3[x"^(c.,Q)]. (2) 



Here a,7 = x,y,z. In the Matsubara formalism, the 



ari component of the susceptibility is given by 
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where (...) denotes the thermal average 



(3) 



(4) 



As a consequence of the SU(2) symmetry of the Heisen- 
berg interaction, all off-diagonal elements of the suscepti- 
bility tensor are zero and all diagonal elements are identi- 
cal. It is therefore sufficient to calculate x^^(w,Q). The 
trace in (|3]) is taken over a basis of states, and Z rep- 
resents the partition function. Using translational in- 
variance, writing the time evolution of spin operators as 
S^ (r) , and inserting a complete set of simultaneous eigen- 
states of the Hamiltonian and the momentum operator 
into the formula for the susceptibility ([3]) gives 

X^^(c.,Q) = -| / dre-"-^^e-*«" ('-'') 
Jo J 1, 
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^ p-T{e„-e^) i{p„-prr,){l-l') Jif 



(5) 



The sum runs over a complete set of states \n) with well 
defined momentum p„ and energy e„ . The expression for 

Mn,m is 
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FIG. 2. The interband transition for {Q\\,Qi_) — (vr, 7r/2) 
found by exact diagonalization of a Jy = 0.25Jx ladder sys- 
tem of L = 8 rungs. The thermal broadening is much greater 
than 77 = 0.01. 



with the result in the thermodynamic limit, we introduce 
a sufficiently large value for the Lorentzian width i] in ([7]) 
to obtain a smooth function. To observe thermal broad- 
ening of the lincshape, the temperature has to be large 
enough for thermal effects to dominate over the artificial 
broadening due to rj. In Fig. [2l we show some typical 
results obtained by this method at intermediate temper- 
atures T > J±. In section I VIII we compare the results of 
the low-temperature expansion developed in the follow- 
ing to the exact numerical answer for L = 8. 



Af„,„=|(n|S'^o|m)| + |(n| ^i^q |m)| 
+ e^Q- {n\ S^o H {m\ Sl^ \n) 
+ e'^Q- {n\ Sl^ \m) {m\ S^, \n) . (6) 

After performing the Fourier transform and analyti- 
cally continuing to real frequencies, Equation ([5]) reads 



x^^(^,Q) 



L 



Y^ 



-/3e„ 



-/3e„ 



2 "^ UJ + in + Cn — Cr^ 
n,m ' 



^Q|l+Pii,P. 



,M„ 



(7) 



II. DIAGONALIZATION OF SHORT CHAINS 

For small systems, we may calculate a basis of simulta- 
neous eigenstates of the Hamiltonian and the momentum 
operator numerically using a standard exact diagonaliza- 
tion (ED) package. This allows the spectral sum in Equa- 
tion d?]) to be evaluated. As a ladder of L rungs has a 
Hilbert space of dimension 4^, this method is only fea- 
sible up to L = 8. The numerically calculated DSF for 
such small finite systems is obtained as a sum over delta 
functions in frequency. In order to facilitate comparisons 



III. LOW TEMPERATURE EXPANSION 



In what follows, we use the fact that for J|| ^ J± 
states can still be labelled by their triplon number for 
J|| = 0, although it ceases to be a good quantum number 
for J|| ^ 0. Subsequently, we will refer to the pertur- 
bative eigenstates as "r-particle states" \^r), where the 
terminology indicates that they reduce to r-triplon states 
when J|| is taken to zero. Here, jr is a complete set 
of quantum numbers uniquely identifying the state un- 
der consideration. Using this notation, we rewrite Equa- 
tion ([7]) as 
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Er,s - 7 E 



3-0e, 
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2 ^ — ' oj + tr] + e^^ — e-y^ 



-^Qll+P-lr.,P',s^'rr,Js 
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For sufficiently small J|| <iC Jj_ , we may associate a formal 
temperature dependence with Er,s and Fj.^s 
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Equation ^ becomes 

^7^,7. = 2 1(7,1 Sl^ |7.)|' (1 + (-l)'^-^cos(gi)) , (10) 
because due to the leg exchange symmetry 

(7.1 5o%(0) |7s) = (-1)'--^ (7.1 SUQ) |7s) . (11) 

The quantities Er^s and Fr^s as well as the partition func- 
tion Z diverge in the thermodynamic limit. We therefore 
reorder the spectral sum in the spirit of a linked-cluster 
expansion following Ref. [3j]. To do so, we express the 
partition function as 



2_^Zn, 



(12) 



where Z„ is the contribution of n-particle states. It is 
furthermore convenient to combine quantities with the 
same formal temperature dependence as 



tfr,; 



Er 



F, 



(13) 



We then define the quantities 



Co 



i=i 



Go, 



Ci — Gi 



C2=G^ 



2,0 + 

oo 

E 

m=2 



7 . (Gij - ZiGo,j-l 
J = l 

G2,i — ZiGifi 



( G2,m — ZiGi^m-l + (-2^1 — ■^2)Go,i 



Cs = . . . (14) 

The C„ are the sums of all cluster functions with the 
same formal temperature dependence. Hence we obtain 
by construction that (as the triplon bandwidth is small 
compared to the triplon gap) 

C„ = 0(e~^"-^^). (15) 

We can then re-express the spectral sum in ^ as 



x^'i^.Q) 



-. oo 



Fr. 



r,s— 



= Ec«- 

n=0 



(16) 



Now we postulate that C„ are finite in the thermody- 
namic limit and (J16p constitutes a low-temperature ex- 
pansion. This assumption is valid in the limit of non- 
interacting dimers J\\ = 0. We furthermore verify it by 
explicit calculation for the leading contribution Ci for 
J|[ ^0. By virtue of the existence of a spectral gap A the 
contribution of C„ is seen to be proportional to e~"^/"^, 
so that p6p constitutes a low-temperature expansion in 
the small parameter e~ '"^, which can be viewed as the 
density of triplons in the state of thermal equilibrium. 



A. Divergences 



As we will see, the expansion ([16)) exhibits "infrared" 
divergences at 

A. w — > ±e((5||). These occur in the "interband tran- 
sition" terms Gj.j+i. 

B. w — > ±2J|| sin((5||/2). These occur in the "intra- 
band transition" terms Gj,j. 

In order to deal with these divergences, we need to sum 
up an infinite number of terms in the low-temperature 
exp ansion. This can be done by following Refs. [23 and 



1. "Interband" processes 

The expansion (|16p contains as the leading term the 
T — Q result, which diverges when the external frequency 
Lo approaches the single-triplon dispersion ^{Q\\) like 



1 



(w -f iriY 



\Q\\) 



(17) 



This corresponds to the coherent propagation of a single 
triplon at r = and leads to a contribution propor- 
tional to 5{lo'^ — e^((5||)) in the DSF. On the other hand, 
for any finite temperature we expect this delta-function 
to be broadened. This is a non-perturbative effect and 
cannot be captured in any finite order in the expansion 
(rT6|) . The fact that a broadening occurs emerges from 
the occurrence of "infrared" divergences in P^ . i.e. sin- 
gularities when the external frequency uj approaches the 
single-triplon dispersion e(Q||). For example, we show 
below that the first sub-leading contribution Ci exhibits 
a divergence 



1 



(w + i77)2-e2(Q||) 



(18) 



We expect the higher terms in the expansion to exhibit 
ever stronger divergences of this type, which need to be 
summed up in order to obtain a physically meaningful 
result. This can be achieved by employing a self-energy 
formalis m^^'^^ . To deal specifically with the divergence 
at llT^ = e2((5||), we divide the expansion (|TB| for the 



susceptibility into a singular (for w^ 
regular piece as follows 



e {Q\\)) and a 



X^^Lo, Q) = XsT„g,i(w, Q) + x?4(^, Q). (19) 

We then introduce a self-energy Si(aj,Q) by expressing 
the singular contribution to the dynamical susceptibility 
in the form of 



Xsing,l(^,Q) = 



Go,l(c^,Q) 



l-Go,i(c^,Q)I]i(w,Q) 
Go,i(c^,Q) + G2 (w,Q)I]i(^,Q)-f. 



(20) 



Here Go,i(w, Q) is the singular contribution to the lead- 
ing term Cq in the expansion ([T51) . Matching (P(7| to 
()16p then yields a low-temperature expansion of both 
Xrcg(w,Q) and the self-energy 



Ei(^,Q)=^I](-''\c.,Q), 



(21) 



where the formal temperature dependence of the n'^ con- 
tribution is 



El"^(c.,Q)=o(e-"'5^). 



2. "Intraband" Processes 



(22) 



In the intraband processes G'jj(w,Q) (j = 1,2), we 
encounter singularities of the form 



4J|2sin2(Q||/2)-(w + i7;)^ 



-J-l/2 



(23) 



We can deal with these singularities by employing a self- 
energy formalism in a way completely analogous to the 
way we proceeded for the interband processes. This re- 
sults in a two-self-energy formalism for the dynamical 
susceptibility. We express x^^('^,Q) as a sum of three 
terms 



x'^tuM) 



Asme 



.l(w, Q) -H Xsfng.2(^, Q) + X?cg(^, Q), 

(24) 
where xlL^,i[^M) and Xsing,2(^'Q) denote the con- 
tributions of all terms singular for w^ — >■ e^((5||) and 
w^ — > 4J|? sin^(Q||/2) respectively. The contribution 



vsing 



2(w, Q) defines a self-energy I]2(w, Q) by 



XsTng.2(W:Q) 



Gi,i(a;,Q) 



l-Gia(a;,Q)E2(w,Q) 

Gi,i(w, Q) -f G?.i(c^, Q)E2(c^, Q) + 



(25) 



Matching the expansions ([25|) to the low-temperature ex- 
pansion for Xshig 2('-^' Q) generates a low-temperature ex- 
pansion for the self-energy E2(w, Q). 



IV. EXCITED STATES IN THE LIMIT OF 
WEAK INTERDIMER COUPLING 

A. Single triplon excited states 

We start with the Hamiltonian ([T]) . Ho is the dominant 
part of the Hamiltonian, which describes L uncoupled 
dimers. The eigenstates of "Hq are tensor products of 
singlet and triplet states at sites n — 0, . . . , i — 1. The 
unique ground state of "Hg is thus a series of singlet states 
on every site n. There are 3L degenerate first excited 



states that consist of i — 1 singlets and one triplet. We 
treat "Hi as a perturbation to "Hq and construct a basis 
for one- and two-particle excited states. 

We define an operator dairn), which creates a triplet 
at site a with z-component of spin m when acting on the 
ground state |0). Single particle states with a definite 
value of momentum p that carry spin-1 are constructed 
as 



L-l 



\p,m) 



4= E ^''"d,,im) |0) 



Vl 



(26) 



With periodic boundary conditions Sl = Sq, trans- 
lational invariance makes momentum a good quantum 
number and the above states are orthogonal, which en- 
ables us to use non-degenerate perturbation theory to 
calculate the single particle energy shifts. To first order 
in a = J\\/J±, the dispersion is given by 

Cp = J_L + J||COs(pa||), (27) 

where an is the separation between the dimers. Imposing 
periodic boundary conditions leads to the quantization 
of one-particle momenta 



JpL 



1. 



(28) 



B. Two-triplon excited states 



We now construct a basis of two particle states in which 
Hi is diagonal. To lowest order in a, the two-particle 
states can be written as as 



L-la-l 



\pi,P2,S,m) =Afs(,Pl,P2) ^^i^aAPUP2)(l)aT l^) , 

(29) 



a=l fc=0 



where 



iS.m 



E $fnr:rn.rfa(mi)4(m2). (30) 



Here $-5^™ are Clebsch-Gordan coefficients. The to- 
tal spin takes values S — 0,1,2 and the normalization 
Ms{pitP2) depends on spin and linear momenta in gen- 
eral. The spatial part of the wavefunction is given by 

^lb{puP2) = e'^P^'^+P^") + vl^^^p^e^^Pi^+P^"), (31) 

where the phase-shifts Af^_^ encode triplon-triplon in- 
teractions. The boundary condition i/'l_i fc(Pi,P2) = 

(— l)"^'0fg(pi,P2), where the sign is due to odd S states 
being antisymmetric, leads to non-trivial quantization of 
two-particle momenta 



i-^rK., 



P2 



— e*^^"^ = g-*P2i 



(32) 



These equations require a numerical solution. Since for 
real momenta Ap_^ is a pure phase, we introduce the 
notation 



<5^ 

Pl,P2 



-*lll«.pj 



(33) 



The normalization of two-particle states is given by 



B. Interband matrix elements 



A/s(pi,P2) ^[L{L- 



1) 



L- 



Si n(^(pi-P2) -'^pipj 
sin(i(pi -_P2)) 



-1/2 



(34) 



The two-particle states have degeneracy 3^(2)- A basis 
of the two-particle subspace in which Hi is diagonal is 
constructed by requiring that 

^2^1 |pi,P2,S',to) = (epi +epj \pi,p2,S,m) . (35) 

Here V2 is the projection operator onto two-particle 
states. This leads to a condition on A^^p^. When the 
triplets in the sum (j27p are not on adjacent rungs, this 
condition is satisfied for any A. Considering the case of 
neighbouring triplets, we find 



Pl,P2 



4I 

Pl,P2 



l^p-i(pi+P2) +2e"'P2 



Ai 



l + e^' 


(P1+P2) 


+ 2e 


-ipi 


1 + e-' 


(P1+P2) 


+ e- 


-ip2 


l + e-' 


{P1+P2) 


+ e- 


-ipi 


1 + e-' 


{P1+P2) 


— e~ 


-ip2 



1 + e^*(Pi+P2) _ e^*Pi 



(36) 



The procedure for solving Equation p6l) is outlined in 
Appendix [Bj Without affecting the result in the thermo- 
dynamic limit, we simplify the calculation by considering 
L to be even. 



V. MATRIX ELEMENTS 



A. Selection rules 

At low temperatures T ^ J±, the leading terms 
in the expansion ([T6| involve states with at most two 
triplons (in the aforementioned sense that the corre- 
sponding states reduce to states with at most two triplets 
in the J|| = limit). In the following we compute the ma- 
trix elements which link 0, 1 and 2-particle states. 

The operator Sji acts on a single site, thus changing 
the triplon number by An = or 1. To first order in a, 
Jii mixes states with those differing in triplon number by 
An — ±2. As however we will only consider the modulus 
squared of the matrix elements, this correction is only 
relevant in the case that the matrix element is non-zero 
to leading order. The rule remains valid. S^ i conserves 
the total S^ which we have used to label states, so AS^ = 
0. The total spin S has to obey the triangle rule. The 
operator under consideration is a vector, thus \AS\ < 1 
and a transition where S' = in both initial and final 
state is forbidden. As the operator is acting on a single 
site, when AS* = the S^ — states have a zero matrix 
element. 



The matrix elements will be expressed in terms of 
Us{p,Pi,P2)- There are several cases to consider for each 
of the types of solution listed in Appendix [B] and their 
respective contributions are shown in Appendix IC II The 
form for a real solution is 

Us{p,Pi,P2) = LNs{pi,P2)e~^^n.v2e^'is 

sin(i(p-Pi + '^pi,P2 -'^S)) 
sin(i(p-pi)) 

sin(i(p-p2-<5pi,p, -T^S)) 



3in(i(p-p2)) 



(37) 



We also calculate the perturbative correction to the ma- 
trix elements to order O (a) in Appendix IC31 The rele- 
vant matrix elements are given in Table HI 



C. Intraband matrix elements 

In the two triplon sector, transitions are possible be- 
tween most combinations of states listed in Appendix |B] 
The full list is shown in Appendix IC2I The result for 
transitions between real states is 



Ws',s{p'i,p'2.Pi,P2) = 

L^US{PI,P2)MS'{P'1,P'2) 



l(« 



2 ^ PI ,P2 p\ .P 



e.„'+(s'-s)^) 



sin(i(pi - K - 5l.,p2 + S^.^,,,^ - {S' - S)n)) 





sin(i(pi-pi)) 






sin(i(pi - 


-P2-SLp2~^pip',- 


is'- 


-S)n)) 




sin(i(pi-p^)) 






sin(i(p2 - 


-Pi+^Lp2+^pip^,- 


(5'- 


-S)7r)) 




sin(i(p2-pi)) 






sin(i(p2 - 


-P2+'5f„P2-Cpi- 


is'- 


-S)n))l 




sin(i(p2-P2)) 







(38) 



In the cases that either of the momenta in the first state 
equals either of those in the second, the corresponding 
fraction needs to be replaced by 

-(L-l)e^^^'--^<i-^-^^'-^^"\ 



VI. SPECTRAL REPRESENTATION AND 
RESUMMATION 

The leading contributions to the low-temperature ex- 
pansion for the dynamical susceptibility are given by 
Go,i = Eq^i + Fi^o- Using the matrix elements from Ta- 



TABLE I. Non-zero matrix elements of S'|o to order a. For the definitions see (|37p . p8|) and (|C16|) . 



(0|SJ,o|p,0) 


(-l)''''i7L(l-fc°«(P)) 




(p',±l|S;,o|p,±l> 
(pi,P2,0,OlS|,o|p,0> 


±2^ 


)-fV^,(p,Pl,P2)) 


(-1)' + Vt2X^ (t^0(p,Pl,P2 


(pi,P2,l,±ljS|,o|p,±l> 


±(-l)'\/5^(t^l(P'Pl'P2) 


-f1/l(p,Pl,P2)) 


(pi,P2,2,0|S|,o|p,0> 


(-l)'VS(t^2(P.Pl'P2)- 


fl^2(p,Pl,P2)) 


(pi,P2,2, ±l|5|,o|p,±l> 


(-l)'V5^(t^2(P'Pl'P2)- 


fV2(p,pl,P2)) 


(pi,p^,l,±l|S|,o|pl,P2,l,±l> 


T4r2-W^l,l(Pl,P2,Pl,P2) 




(pi,P^,2,±2|S|,o|pi,P2,2,±2) 


T2l^W^2,2(pi,P2,Pl,P2) 




(p'i,p^,2,±l|S|,o|pi,P2,2,±l) 


±4l2-W^2,2(pi,P2,Pl,P2) 




(p'i,P^,2,±l|S|,o|pi,P2,l,±l> 


±4i7W'2,l(pi,P2,Pl,P2) 




(p'i,P^,2,0|S|,o|pi,P2,l,0) 


2^^2^2,l(Pl,P2,Pl,P2) 




(p'i,P^,l,OlS|,o|pi,P2,0,0) 


2r^\/ 3^1.o(Pl'P2,Pl,P2) 





blelH we find that to order a we have 
(1 - cosQ^) 



Go,: 



(1 — acosQ||) 
1 1 



w + i7/-eQ w + iry + eg 



(39) 



These give rise to a delta function peak located at the 
one-triplon excitation energy. The intraband term G\,\ 
is given by 



Gi,i 



(1 -fcosQj 
2L 



E 



-/3ep 



-/3eQ|| 



w + i7] + ep-eQ||+p 



(40) 



Similarly, we find the interband terms 



Gl,: 



(1 — cosQj 
4L2 



■E(^ 

Pl>P2 



-/3eQ|| 



1 



^w + ir; + eQ||+pi+p2 ^ f^pi ^ '^P2 
1 \ 



uj + m; + epi + ep 



^Q||+Pl+P2 



xE^(|f^5-«C^s^s|). 



(41) 



The sum over pi,p2 is taken over all momenta that 
satisfy the boundary conditions p2p , and these momenta 
depend on 5*. The leading term in 6*1^2 scales with L, but 
cancels against the "disconnected" contribution ZiGo,i- 
The low-temperature expansion of the dynamical suscep- 
tibility now takes the form 



x""(c^,Q)«Co+Ci+C2, 



(42) 



where 



Co(ii^,Q) ~ Go,i, 
Ci(l^,Q)«Gi,o + Gi,i + (Gi,2 

C2(c^,Q)«G2,2-^lGi4. 



ZiGq: 



(43) 



Here Zi — 3 ^ 6"^*^^ is the single particle contribution 
to the partition function. We note that in C2 we only have 
taken into account the intraband processes. We observe 
the following divergences in C„: 



C„(w,(3||) ex 



(c^+»r,)^~^^(Q||) 



1+n 



^1/2 



where we have defined 

e(/c) = 2J||Sin(Q,|/2). 



^2 



(44) 



(45) 



The first (second) kind of singularity is seen to be present 
in C„ for n = 0,1 {n = 1,2). We expect (gl]) to hold 
for n > 2 as well. Following the procedure set out in 
Section lillAl we define 

Xshig,2 ~ Gi,i + (G2,2 — ^lGi,i), 

XsTng.i « Go,i + Gi,o + (Gi,2 - ZiGo,i). (46) 

The leading orders in the low-temperature expansions of 
the self-energies then take the form 

Si(c., Q) - G^^,{uj, Q) [Gi^2(c^, Q) - ^iGo,i(c<., Q)] , 
5]2(w, Q) = Grj(w, Q) [G2,2(c^, Q) - ^iGi,i(w, Q)] . 

(47) 



Our approximate result for the DSF is then 
5*^^(0;, Q) =- limi ^ 



ri'-i'O TT 1 — e^'3" 

Gi,i(^,Q) 



l-Gi,i(w,Q)I]2(c^,Q) 
1-Go,i(c^,Q)Si(l^,Q) 



(48) 



VII. RESULTS AND DISCUSSION 

In order to present expUcit results, we choose a = 0.1 
and perform numerical calculations on a system of L = 
1000 dimers. Doubling the number did not change the 
results significantly. The limit 77 — >■ is approximated by 
choosing a value larger than the spacing of the momen- 
tum values due to finite size, which is of order 0{^J\\), 
but small compared to the thermal broadening Jye"'^'^^, 
so that the shape of the response is not changed signif- 
icantly. One problem we encounter is that to the order 
in J\\/J± we are working in, the bound state contribu- 
tions to Ci give rise to sharp peaks for kinematic reasons. 
These features will be suppressed once higher orders in 
perturbation theory are taken into account, even if we 
do not sum higher order terms in the low-temperature 
expansion (which would lead to a further broadening). 
Given that the sharp bound state peaks are an artifact 
of the order in perturbation theory considered, we choose 
to suppress them in the various plots by specifying a suf- 
ficiently large broadening ij = 0.01. This also facilitates 
comparison to the ED results. The choice of Q± affects 
the mixing between the intraband (oc cos^ hQ^) ^^^ in- 
terband (oc sin hQ±) responses. Hence plots are given 
for Q± — 7r/2, where both types of transition are allowed 
with equal weight. 
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FIG. 3. The interband transition for (Qy , Q±) = (tt, tt/2) and 
L = 1000 sites. The asymmetry grows as T increases. 



resummed 
Lorentzian fit 




FIG. 4. The resummed interband transition lineshape for 
T = 0.4J_i, Q± = 7r/2, Qy = n, rj = 0.01 and L = 1000 
together with a Lorentzian best fit demonstrating the asym- 
metric lineshape. 



A. Broadening of the triplon line 

We first consider the temperature evolution of the 
triplon line. At T = the DSF features a delta func- 
tion line following the triplon dispersion. In Fig. |31 
we plot S'^^(a;,Q) as a function of frequency for wave 
vector Q = (7r,7r/2) and temperatures in the range 
0.2 J± < T < 0.4J^. We see that the line broadens asym- 
metrically in energy as the temperature increases. On the 
other hand, at sufficiently low temperatures we expect 
the lineshape to be well approximated by a LorentziariF^. 
In Fig. 21 we show a comparison of the actual result to a 
Lorentzian fit 



5Lor(w,Q)=A(Q) 



1/T0 



(a.-e(Q||))2 + l/r| 



(49) 



Fig. [S] shows the dependence of the asymmetry on (5||. 
The falloff is slower towards the centre of the dispersion. 




FIG. 5. Dependence of the interband transition at T = 0.5 J_l 
on (5||. Q± is fixed at 7r/2 and L = 1000. The graphs are 
offset by n for Qy — nvr/l. 
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FIG. 6. A comparison of the resummed spectral function for 
T = 0.4 Jx, Q± = tt/2, Q|| = tt, »7 = 0.02 and L = 1000 to 
the ED result for a L = 8 system. 
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FIG. 7. The intraband transition at a series of temperatures. 
Q± = 7r/2, Q,| = TT, r? = 0.01 and L = 1000. 



In order to establish the temperature range in which our 
low-temperature expansion provides accurate results, we 
compare (P5|) to numerical results obtained by a direct 
diagonalization of the Hamiltonian for short chains. To 
obtain a continuous curve for the DSF, we convolve the 
numerical results with a Lorentzian of width rj = 0.02. 
Fig. [S] shows such a comparison for T = QAJ±, Q± = 
7r/2, Q|| ~ TT and L = 1000. We see that there is good 
agreement between the two methods. 



B. Finite temperature resonance at low^ frequencies 

As in the state of thermal equilibrium there is a finite 
density of triplons, incident neutrons can scatter off them 
with energy transfers small compared to the gap. Ac- 
cordingly at finite temperatures there is a spin response 
at energies uj ^ 0. To leading contribution to this "intra- 



band response" is 



l + cos{Q^) e 



-l3(J^u/2) 



TT 1 - e-'3" 
X cosh 



27r 



^e\Q\\)~u^ 



^ /3cot(Q||/2) r— — ^\ „ . 2,n A 



(50) 



where e(Q||) is given by pS)) . This contribution con- 
tains square root singularities for u; — >■ ±e((5||), which 
get smoothened once we resum terms following Sec- 
tion llll A^ In Fig. [7] we plot the DSF at low frequencies 
for several temperatures in the range 0.2 J_l <T< 0.4 J_l. 
We see that the integrated intensity increases with tem- 
perature, while a strong peak at w sa e^Qii) remains. This 
is very similar to what happens in the spin-i Heisenberg- 
Ising chairt^i, where this feature was first predicted by 
Villain^. 



C. Summary 



In this work we have determined the low temperature 
dynamical structure factor of the two-leg spin-i Heisen- 
berg ladder in the limit where the leg coupling is weak 
compared to the rung exchange. We have shown that 
the sharp delta-function line following the triplon disper- 
sion at T = gets broadened in an asymmetric way at 
T > 0. The dominant processes at low T involve scatter- 
ing from one-triplon to two-triplon states in the presence 
of a "thermal background" , as described in Section IIIII 
We have also determined the temperature activated con- 
tribution to the DSF at low frequencies. Here the dom- 
inant processes at low T involve scattering between dif- 
ferent two-triplon states in the presence of a "thermal 
background" . Our analysis is based on the method de- 
veloped in Ref. [2J| for the case of the alternating spin-i 
Heisenberg chain. We have gone beyond Ref. '24l| in two 
important aspects. Firstly, we have taken into account 
all perturbative corrections to the various matrix ele- 
ments to order O {J\\/J±)- This establishes that higher 
order perturbation theory in J\\/J± can be combined 
with the low-temperature expansion of Ref. [2J]. Sec- 
ondly, we have included the order O (e"^^''^) correction 
6*2,2 ~ ZiGi^i to the intraband contribution. This allows 
us to describe the low-frequency temperature induced 
"resonance" in a significantly larger temperature window 
and demonstrates the difficulties encountered when deal- 
ing with higher orders in the low-temperature expansion. 
It would be interesting to compare our results to exper- 
iments on ladder materials. Perhaps the best candidate 
is (C5Hi2N)2CuBr4, which is a highly one-dimensional 
two-leg ladder material with a ~ 0.256^"—. Experimen- 
tal studies of the temperature evolution of the DSF for 
this material are under way^. 
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Appendix A: Linked-cluster expansion for Ji 







For Jji = we are deahng with an ensemble of uncou- 
pled dimers. The dynamical susceptibility can then be 
calcuated by elementary means in the Matsubara formal- 
ism. After analytic continuation we obtain 

-I3J± 



X^^(^>0,Q) 



J± 1 



1 — C0s((5j 



(Al) 



2 1 + ie-f^J^ [lo + i0)2 - J2 

The temperature dependent factor can be expanded at 
low temperatures 

1 _ p-PJ± 

^, = 1 - Ae~^'^^ + 126-2^-^^ H- . . . . (A2) 

We have calculated the first few terms of the low- 
temperature expansion (J16l) by working in a product ba- 
sis of dimer triplet and singlet states. The leading con- 
tribution is 

J± 1 - cos{Q±) 



Co — Go.i — 



(A3) 



2 (w + i0)2-j2' 

which correctly reproduces the T = limit of (|Aip . The 
next term is 

Ci = Gi.o + (Gi,2-^iGo.i). (A4) 

We find by explicit calculation that 

Gi,2 - ^iGo,i - -Se^^'-'^Go,!. (A5) 

This results in 

Ci=-4e-^^^Go,i, (A6) 

which correctly reproduces the first subleading term in 
(IA1[) . The next term is 

C2 =(G2,i — ZiGifi) 

+ (G2.3 - ^iGi.2 + {Zl - Z2)Go.i). (A7) 

We find that 

G24 — ^iGi_o ~ — 3e ^Go,i, 

G2,3 - ^iGi,2 + {Zf - Z2)Go,i - 9e-2'3''^Go,i, (A8) 

which gives 



Co = 12e 



-HJi 



Go,i. 



(A9) 



This correctly reproduces the second subleading term 
in (jAl|) . We note that in the limit J|| = the low- 
temperature expansion (J16p is well defined and does not 
suffer from the kind of "infrared" divergences present for 
J|| ^ 0. This is as expected since the spectral function 
of the full result (jAl|) features a sharp delta-function line 
even at T > 0. 



Appendix B: Solutions of the BAE 

1. Real solutions 

To find the two-triplon momenta allowed by the quan- 
tization condition p2l) . we follow the approach outlined 
by James et ali^. We choose a suitable branch cut such 
that the solutions are enumerated by 



Lpu2 = TiH-Ap 



2n 



h. 



l + (-l)' 



, (Bl) 



where /i^2 are integers used to parametrize the equation. 
This gives L{L — l)/2 possible solutions. To satisfy pi > 
P2 we need /i > l2- In the case of /i ^ I2 this is easily 
solved numerically, although care must be taken not to 
double-count solutions. 

One must be careful with those solutions where the 
phase shift is zero. The momenta are then equal to the 
single triplon momenta and so the matrix elements can 
be of order 0{L). These solutions occur only for S ~ 
or 2. For these states the normalization is 



Ms - [L{L - 2)] 



(B2) 



The procedure above does not identify all real solu- 
tions in the S = sector. The remaining roots are found 
following Ref. [401 ■ ^'^^ large systems Equation (|B1[) has 
solutions where /i = l2- Whereas the trivial solution 
Pi — P2 is forbidden by the Pauli principle, another so- 
lution appears very close to the trivial one. Due to the 
proximity of the two zeros, the numerical solution of the 
equation is difficult. One method is to rewrite the BAE 
as a single equation in x — pi — p2 and to then divide by 
X to eliminate the trivial zero, after which the root finder 
converges reliably on the desired solution. 



2. Bound states 

There also exist complex solutions pi^2 = x±iy, where 
the amplitude decays exponentially as a function of the 
separation of triplons, corresponding to bound states. 
For these the S-matrix elements are real, and equation 
becomes 



gi^L^-yL 



+ (-1) 



g2cos(x)-H(2-|(5'+l))e-^ 
2cos{x) + {2-§{S+l))ey 



= 0. 



(B3) 



For each x = mr/L there may exist a zero, and the num- 
ber of solutions scales as L. The matrix elements for 
these roots require special treatment and are given as 
previously in terms of 






+ L- 



2sinh(y) 



(B4) 



10 



3. Singular solutions (type I) 

At this point we still miss 4 solutions, which occur at 
singularities of the quantization conditions. Such a solu- 
tion was described for the spin-i XXX model in Ref. pOj . 
For each S sector there is a solution at pi^2 = 7r/2 ± ioo, 
corresponding to a vanishing S-matrix eigenvalue. By 
introducing a twist angle the quantization conditions 
become 



At _e**/2 = i~lfe'^P\ 



gi<t> _ g«i(pi+P2) 



Pl,P2'- 



(B5) 



This renders the momenta finite, but they cease to be 
complex conjugate to one another. Normalizing the wave 
function and then taking the limit — >■ we obtain 



V-f.h = {-l)HSa-l,b - hlfSa^L-lSbM). 



(B6) 



It can be verfied by direct calculation that this gives an 
eigenstate of the Hamiltonian. The normalization of the 
state is 

J\fs^L-i. (B7) 



4. Singular solutions (type II) 

Finally, there is another singular solution in the 5 = 
sector with pi — P2 — tt. This solution gives rise to an 
eigenstate despite the fact that the two momenta are the 
same because the phase shift is ill-defined. Again the 
limiting wave function can be calculated by introducing 
a twist angle, normalizing the state and then taking the 
twist angle to zero. The result for the wavefunction and 
its normalization is 



^% = (-1) 



a+b 



AAn 



L{L-1) 



(B8) 
(B9) 



Appendix C: Matrix elements 

1. Interband matrix elements 

The interband matrix elements for the different types 
of solution are as follows: 

A. Real solutions with zero phase shift: 
Usip,Pi,P2) = -L^s (L idp,,p + Sp,,p) - 2) . (CI) 



B. Bound states: 
Us{p,pi,pl) =LAfs{pi,pl)e''^ 



cosh(y) — cos(a; — p) 

nS ^ 



i^ + ^P^.Pt)cosix^p-—) 



C. Singular solutions (type I): 

Usip) = 2iLAfssinip). 

D. Singular solution (type II): 

Uoip) = LMo. 

2. Intraband matrix elements 



(C2) 

(C3) 
(C4) 



The intraband matrix elements are as follows for tran- 
sitions between different types of states in the two triplon 
sector: 

A. Real — > Bound: 

Ws'.,s{p'i,p'*i,Pi,P2) = 

i'AAs(pi,P2)AA5'(K,p'I)e*(^'~^^"e^^'i-2 

(<.P'i^" + ^"") '^"^(('^ - ^')f + ¥Lp.) 

cos(pi — x) — cosh(y) 
(<,p-» + l)cos{x-pi jS-S')^ ^ l^g^,,J 
cos(pi — x) — cosh(y) 

cos(p2 ~ x) ~ cosh(y) 
(Ag^^,. + 1) cos(a; - p2 - (5 - 5')f + ^^LpJ 



+ 



cos(p2 — x) — cosh(y) 



B. Real — >■ Singular (type I): 

Ws\S[^,^,Pl,P2) = 

L2AAs(pi,P2)A/'s'e*^'--+(^'-^)"^ 



(C5) 



X 2 



cosL-^-(5'-5)| 



cos P2 + 






(C6) 



C. Real — >■ Singular (type 11): 

Wo,i(7r,7r,pi,p2) = 

- iLWi(pi,P2)A/'oe3''pi.P2 



Pl,P2 



tan h tan — 

2 2 



(C7) 
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D. Bound — > Bound: 



a. Ground state corrections 



Ws'.,s{p'i,p'l,Pl,P*i) = 

L^Ms{pi,plWs'{p[,p'l)e'^''-'''^ 

1 — e'^^'-^^^'^^y^y' 

AS' _ (_-,\S+S' AS p-i{x-x')+y-y' 

-y+y' P'l.P'l '^ '' ^PuPr 

I _ gi(x-x')-y+y' 

aS _ (_^\S+S' aS' -i{x-x')-y+y' 

j-y' ^Pl,Pl y ^' ^P'l.p'l^ 

I _ gi{x-x')+y-y' 

AS AS' _ (_-,\S+S' -i{x-x')-y-y' ■ 

y'^PuPiP'i.p'l y '-> ^ 



aV+y 



\ _ pi(x-x')+y+y' 



(C8) 



E. Bound — > Singular (type I): 



Ws',s{i^,l^,Pi,P2) 






(C9) 



F. Bound — > Singular (type II): 



Wifi{pi,pl,TT,Ti) = iL^AfiAfo 



sin(a;)(l + Ai^p.) 



cos(x) + cosh(y) 



(CIO) 



G. Singular (type I) — > Singular (type I): 



Ws'A\. \. f , \) ^ '2L^AfsAfs'5s,S'. (Cll) 



H. Singular (type I) — > Singular (type II): 



WoAt^,t^, -, -) = -2LWiAAo. (C12) 



3. Corrections to the interband matrix elements to 
first order in a 

We expand the states to first order in a and calcu- 
late the corrections to the matrix elements. Firstly we 
note that "Hi can only induce transitions between states 
where the particle number differs by at most 2, since 
each term in the sum only acts on a pair of adjacent 
rungs. Secondly, the Hamiltonian is symmetric under 
leg-exchange, while a state \js) with s particles picks up 
a sign of (—1)^^". This implies (7^1^1173) = if |r — s| 
is odd. Hence the only contributions to first order are 
from states with a particle number different by 2. 



These are given by 



72> 

= 10) + ^a 



~2J 



L-l 



J2<^'aka\0)+O{a^). (C13) 



a=0 



b. Single particle state corrections 

There are the following contributions from three- 
particle states 



(1) _ \r- (73 1 "Hi \p,m) 



|p,m)*^^ =^ 



|73> 



a^i 



-2J, 



L-l 



173; 



^ a— hi^a.a—1 



(C14) 



c. Two-particle state corrections 



In the two-particle sector there are contributions from 
four-particle states and from the ground state. The for- 
mer do not contribute to any of the matrix elements used 
in the subsequent calculation to first order, so they are 
not calculated. As the Hamiltonian conserves S and to, 
only the \pi^P2t S — 0,m = 0) state will have a correction 
from the ground state. For real solutions this is given by 



bi,P2,o,o)(^' 



(0|Hibi,P2,0,0) 
2J, 



|0) 



"Pl+P2,0 



L V3 
L-l 4 



a(e*Pi - 1) 



(C15) 



Hence the corrections to the matrix elements are: 
A. Real solutions: 
Vs{p,pi,P2) = LUs{pi,P2)e'^^^^-''^e^^'^ 



sm{^{p-pi+S^^p^ -t:S)) 
sin(i(p-pi)) 

sin {}{p-p2- (5^1 ,p2 - ttS)) 



pi + 2p2-p 



P1.P2 ■•-'J J I 2pi+p2-p 

'^-^ cos 



sin(i(p-p2)) 
3o5o 2 cos cos 



2 

Pi +P2 



-LSp^+p^fi cos 



2pi 



6" 

Pl,P2 



(C16) 



B. Real solutions with zero phase shift: 



2cos(— ^) - L6p,^p\ cos(^ P2) 



2cos( )-Ldp^^p ) cos(^ pi) 



+ 5so& cos( — - — ) cos( — - — , 



C. Bound states: 



(C17) 



D. Singular solutions (type I): 



Vo,2(p, I' |) = iLNs sin(2p). 



E. Singular solutions (type II): 
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(C19) 



Vs{p,Pi,pI) = ■^LJ\fsipi,pl) 



HSs,oS2x,oL + 1) cosx{e~y + A^,,p.e«) 



+ 



+ 



cosh(y) — cos(x — p) 
{e-y + A^^p.ey){cosip - ttS) - cos{p - 2x)) 

cosh(y) — cos(a; — p) 
ie~^y + AS^^^,e^v)cos{x-7rS) 
cosh(2/) — cos(a; — p) 



(C18) 



Vo{p,n,Tr) = LNs 



2 sin 



2P 



(C20) 



The matrix elements with their corrections can be found 
in Table HI There is also a non-zero correction to the 
matrix elements contributing to £'0,2 and Fo_2, but since 
this term is zero to leading order the correction to the 
modulus squared of the matrix element is only second 
order. 
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